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X ■ ABSTRACT 

We formulate the particle conjugation operation and its convenient realization as G- 
parity in the framework of several chiral soliton models. The Skyrme model, the Skyrme 
model with vector mesons and the chiral quark model are specifically treated. The vector 
and axial vector currents are classified according to their behavior under G-parity. In the 
soliton sector particle conjugation constrains a priori ambiguous orderings of operators in the 
space of the collective coordinates. In the Skyrme model with vector mesons and in a local 
chiral model with an explicit valence quark this classification scheme provides consistency 
conditions for the ordering of the collective operators appearing in the 1/Nc corrections to the 
nucleon axial charge and the isovector magnetic moment. These consistency conditions cause 
the corrections obtained from an ordinary perturbation expansion to vanish in the context 
of the collective quantization of the static soliton configuration. This conclusion presumably 
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applies to all local effective chiral models. 
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1. Introduction 



In this paper we shall discuss some theoretical constraints on the vector and axial vector 
current matrix elements of the nucleon predicted in various kinds of chiral soliton models (for 
review articles see |^). One might, at first, think that this is a subject which has been 
completely exhausted. However, there are, at least, two recent problems: 

i) the "too small qa problem" 

ii) the "proton spin puzzle" ^ 

where such considerations appear crucial. These problems involve current matrix elements 
which are of subleading order in an appropriate perturbation approach to the collective 
Hamiltonian describing the nucleon. As such, they turn out to involve operator ordering 
ambiguities. Symmetries like iso-spin invariance do not seem to be able to offer guidance 
on this question. Here we point out that particle conjugation symmetry it does in 

ordinary field theory, provide operator ordering restrictions. For example, charge conjugation 
invariance requires the QED interaction to have the ordered form ^eA^ (^'^^^> — \l/-^7j\l/"^^ . 

Recently, the very interesting observation has been made 0-0 that the too small value 
of the neutron beta decay constant qa predicted in many chiral soliton models might be 
dramatically improved by including 1/Nc corrections. This seems satisfying both because it 
agrees with the quark model result in which qa should behave like {Nq + 2)/3 so sizable 
corrections are expected and because the corrections fit the pattern expected from unitarity 
constraints on pion-nucleon matrix elements |TT|] (which imply that the structure of the 
leading and next-to leading order terms is identical). In spite of these pleasing features, we 
will see that the imposition of proper ordering constraints requires that these 1/Nc corrections 
vanish in certain chiral soliton models. 

Actually, a different problem with these corrections has already been noted |^. The 



equation of motion for the pion fields, which yields the PCAC (Partially Conserved Axial 
Current) relation, is derived from the static energy functional of the soliton. This functional 
is of the order Nc without any subleading corrections. Thus this equation of motion cannot 
account for 1/Nc corrections in the divergence of the axial current. The pion field develops 
a Yukawa-type tail whose amplitude is, due to PCAC, proportional to qa- This allows one 
to not only compute qa by direct integration of the axial current but also from the large 



distance behavior of the pion |13|. The latter, however, is completely determined from the 



classical equation of motion and hence does not contain the l/Nc corrections. Therefore it 



has been concluded that these corrections violate PCAC by about 30% |]T2[. This problem 
was circumvented by elevating the derivative of the axial current including the corrections 
to the equation of motion via PCAC [|1^, 0. However, this treatment is not completely 
satisfactory since it has not proven possible to derive the so constructed equation of motion 
from an action principle. 

We shall, for simplicity, work in a theory with just the two lightest quark flavors present. 
The relevant currents at the quark level are the vector and axial vector ones of both isoscalar 
and and isovector type 

K = lill.l^q-l'llllf)^ A; = ^{q^,^,T^q-q^jj:jlT^f), (1.1) 

where q represents the quark field. Notice that all of these operators are hermitean. In 
dealing with the charge conjugation of objects with definite isospin transformation properties. 
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it is often useful to define the concept of G-parity [0. Technically the G-parity reflection 
is defined as charge conjugation followed by a rotation of angle vr around the y-axis in 
isospace. This corresponds to a unitary operator Q = e^'^^^C. Since the G-parity of the pion 
is minus one, the G-parity classification provides a tool to determine whether the number 
of pions emerging from a decay process is even (G-parity= 1) or odd (G-parity= —1). By 
construction the G-parity is a conserved quantity for interactions which are invariant under 
charge conjugation and isospin rotations. For the current operators in ( p..l[ ) we have 

g (v;°, v;, Al A») g-' = (-v;°, v;, aJ, -a;) . (1.2) 

In chiral models the pions are incorporated via the non-linear realization U = exp{i7z ■ 
t/Ztt) where = 93MeV denotes the weak pion decay constant. Since tt — > — tt under the 
G-transformation we simply have 

U{r,t) ^ U\r,t). (1.3) 

One should bear in mind that in all quark soliton models the quark fields are functionals of 
the chiral field so their behavior under conjugation is determined from that of U. 

This paper is organized as follows. In section 2 the behavior of the symmetry currents in 
the simple Skyrme model under G-parity reflection will be studied. This will be generalized 
to models with vector mesons in section 3. In section 4 we will first discuss the particle 
conjugation in the framework of quark soliton models and then apply these results to the 
l/Nc corrections to and /xy. 

Actually the problem of ordering ambiguities in soliton models was already discussed by 
Tomboulis |]14| in the case of a two dimensional field theory. Also in that reference symmetry 
arguments (Poincare invariance) have been employed to resolve operator ordering ambiguities. 
However, we will mainly be concerned with form factors at zero momentum transfer and hence 
do not introduce collective coordinates for the boost of the soliton. Rather we are interested 
in ordering ambiguities for the collective coordinates associated with internal symmetries^. 

2. G-Parity Constraints on the Currents in the Skyrme Model 

In order to introduce the relevant background and notation we start with the simple 



Skyrme model, defined by the Lagrangian [|T^ |T3 

7T 

4^' V"'^ / ' 32e 



£ = ^tr (d^Ud'^U^) + -^tT{[d,U, dM^] [d'^U, d'U^]) . (2.1 



where e is the only free parameter. 

The static soliton configuration of the Skyrme model is given by the celebrated hedgehog 
ansatz 

Uo{r) = exp[iT ■ rF{r)] (2.2) 

which introduces the chiral angle F{r). Substituting this ansatz into the Lagrangian (|2.1| ) 
yields minus the classical energy functional E^. The soliton configuration is then obtained 
by minimizing this functional. In the Skyrme model the baryon number current is identified 
Q with the topological current which classifies the mappings of the type (p^) with the 



"For a two dimensional model operator orderings were considered in ref. with regard to internal 

symmetries. 
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boundary condition F(r) — >0. The baryon number associated with a given chiral angle, 
F{r), is obtained as the spatial integral over the time component of this current 



2 

B = - 



/ drF\r)?>in^F{r) = F(0). (2.3) 

JO vr 



Integer baryon numbers thus correspond to the boundary conditions F{0) = nn with n being 
a positive or negative integer. In this context it should be remarked that in topological soliton 
models Ed diverges as one deviates from this boundary condition. Obviously, one crosses 
from the baryon to the anti-baryon sector by reversing the sign of the chiral angle, F{r), 
since not only the baryon number but also the baryon number density acquires the opposite 
sign. 

In order to project the soliton onto states with good spin and isopin one introduces 
time dependent collective coordinates parametrizing the corresponding rotations. Due to the 
symmetry of the hedgehog ansatz the rotations in coordinate and iso-space are equivalent. 
Hence one approximates the time dependent solutions by [|l^ 



U{r,t)=A{t)Uo{r)A\t). (2.4) 

The collective coordinates are contained in the SU(2) matrix A{t). From the defining relation 
( |1.3| ) we observe that the G-operation corresponds to the reversal of the sign of the chiral 
angle F — * —F while keeping the collective coordinates, A{t), unaltered. According to ( |2.3| ) 
this indicates that in the soliton sector the G-parity transformation actually corresponds to 
a particle conjugation. 

A more transparent parametrization of the collective coordinates and their time depen- 
dence is given in terms of the adjoint representation matrix. Dab, and the angular velocity, 

n, 

Dab = - ti (raAnA^) and A^-A='-T-n. (2.5) 



The collective Lagrangian associated with the isorotating hedgehog ( p.4|) 

L = -Ed + ^a^n' (2.6) 

defines the moment of inertia, a^, which again is a functional of the chiral angle, -F(r)0. Here 
we should stress that the classical mass as well as the moment of inertia do not change when 
the sign of the chiral angle is reversed, i. e. when one goes from the baryon to the anti-baryon. 
Furthermore both quantities are of the order Nc- 

In the process of quantizing the collective coordinates, the spin operator is identified as 
the quantity canonical to the angular velocity ft 

J = § = a.^a (2.7) 

For what comes later it important to note that the spin operators act as right generators in 
the collective space 

[J, A] = ^At =^ [J„ Dab] = ieibcDac. (2.8) 



''For explicit expressions of Ed and in various Skyrme type models the reader should consult ref. 
and references therein. 
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Due to the invariance of the static sohton under combined rotations in coordinate and iso- 
space the isospin is related to the spin via the rotation 

la = -DabJb- (2.9) 

This restricts the possible eigenstates of the collective Hamiltonian 

H^E,,^£,^E.,,'I^±11 (2.10) 

to those with total isospin I = J. Since the isospin operators correspond to the left generators 
one may replace 

Dab -paJb (2.11) 

in the subspace / = J= l/2asa consequence of the Wigner-Eckart theorem. We therefore 
do not need any further specification of the nucleon or anti nucleon wave-functions in the 
collective space. 

Various static nucleon properties are obtained from the matrix elements of the vector and 
axial vector currents (|1 . 1| ) . These correspond in the Skyrme model to the symmetry currents 
which couple linearly to external gauge fields introduced in eqn ( p.l| ). The currents are then 
given in terms of the chiral field and its derivatives. Next the rotating hedgehog configuration 
( |2.4D is substituted and the angular velocity is eliminated with the quantization prescription 
( [2.71 ). As the moment of inertia, is of the order Nq it is obvious that each such substitution 
reduces the order of Nc by one unit. Finally the currents are obtained as combinations of 
collective operators whose coefficients are radial functions which depend on the chiral angle, 
F[r). These combinations are displayed schematically in table |2?T| . Note that due to the static 
nature of the soliton configuration the time (/i = 0) and spatial (/i = i) components of the 



currents behave differently. Actually the form of currents displayed in table ^]T] is the most 
general one at leading order in 1/Nc once proper account is taken the rotational symmetries 
in coordinate- and isospace as well as parity. Additionally one might have a time component 
for the axial singlet current Aq = y4o(r)f • J. However, AQ{r) = in Skyrme type models. 
We have defined the isospin part of the generator of the isoscalar currents to be 1/Nc- The 
proper normalization of the electromagnetic charges is guaranteed by / d^rVi{r) = 1 and 
/ d^rV^i^r) = o? . The explicit forms of the radial functions Vxif)^ . . ^ A^ir) may readily be 
extracted from the review articles ilH. 



We have also listed in table |2.1| , the possible sign change of the currents when one replaces 
F —>■ —F in the explicit forms of the radial functions. According to our discussion above this 
sign change should represent the G-parity of the appropriate current. In every case, this is 



seen to agree with ( |1.2| ) in the underlying QCD theory. 

It is important to note that in the transition from the classical objects Dab and O to 
the quantum operators Dab and J there can be ordering ambiguities in cases where the 
product of two operators is required. Such products appear for the time components of the 
isovector- vector current, Vq, and the isovector-axial vector current, Aq. For the former there 
is no ambiguity because [Dab, Jb] = 0. For the latter we must choose a symmetric ordering 
to preserve the hermitean nature of the current. Using the substitution ( p.llD shows that 



matrix elements of Aq vanish between nucleon states. This has the desired consequence that 
the G-parity violating form factor Gxiq^) in the Lorentz covariant decomposition of the 
matrix elements 



{N'\Al\N)=u{p'\ 



GA{q + + 2M ^^"^^ 



■J- 
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Table 2.1: The structure of the symmetry currents in chiral sohton models, the power of the 
leading term in an 1/Nc expansion as well as the behavior of the radial functions under sign 
reversal of the chiral angle and their G-parity. 



Current 


Structure 


0{Nc) 


F ^ -F 


G-parity 




Vi(r) 









' I 


V2(r)eyfcrjJfc/a^ 


-1 






ya 


-V^{r)D^^J^Ic? = V3ir)IJa^ 





+ 


+ 


vt 




1 


+ 


+ 











+ 


4° 


{Ai{r)6ij+A2{r)firj)Jj/a^ 


-1 


+ 


+ 




A3{r)eijkri {Daj,Jk]/oP' 











{Ai{r)5ij + A^{r)rirj) Daj 


1 







is zero in the Skyrme model (c/. appendix C of ref. |T^). If we had incorrectly chosen a 
different ordering for the operators, Aq would not have vanished between nucleon states and 
a non zero Gxiq^) form factor would have resulted, violating G-parity. This provides an 
illustration of the connection between proper operator ordering and G-parity invariance. At 
this stage we see that a symmetric, or hermitean, ordering of the collective operator is in 
accordance with the G-parity symmetry. 

Returning to the discussion of table ^Tl] we again stress that when going from the baryon 
to the anti-baryon sector, i.e. F ^ —F, the symmetry currents acquire a sign according 
to their G-parity quantum numbers. This just reflects the fact that the number of pions 
(as measured by the power in which F appears) is either even or odd for a given current. 
Of course, all extensions of the model should maintain this property. Furthermore, operator 
ordering ambiguities of a current, which may occur at higher order in the I/Nq expansion, 
have to be resolved in such a way that the phase acquired under F ^ —F is determined by 
the G-parity quantum number of the current. 

This prescription has a fruitful application for the axial-singlet current, which is especially 
interesting in the context of the proton spin puzzle. Not long ago it was proposed that an 
axial-singlet current of the form 

Al = ic tr (udM^Udf^WUd^W) = ic tr {puPf^p'') , (2.13) 

(resulting from an additional term in (p.l|) ) would lead to a non- vanishing matrix element 
between proton states |jl9|. The anti-Hermitian pseudovector is given by = d^C,C,^+^'^dfj^^ 



with ^ denoting a root of the chiral field, i.e. U = The reversal F —F corresponds to 
^ and thus —p^. It is then obvious [^] that, unless the current ( |2.13|) vanishes 



identically, it contradicts the G-parity symmetry, cf. table p.l| . Let us see, however, how an 
improper ordering can lead to a non-vanishing result. As = p^r^ is a vector in iso-space, 
( |2.13|) can be expressed as 2ieabcPlP^fjP'^'^ which vanishes as long as the are considered 



to be classical objects. Under the canonical quantization (|2.7| ) the angular velocity, O, 
which is contained in po, is replaced by J joP'. Then one is tempted to replace eahS^h^c by 
(■abcJbJc/ {oP'Y = iJajip^Y- This in turn would yield a non-vanishing axial-singlet current 

4° = r^F'{r)sm^F{r)r ■ Jr, = -— -rVi(r)f ■ Jf,. (2.14) 
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However, ( p.l4| ) is suspicious on a priori grounds since it changes sign on F — > —F, in dis- 
agreement with the requirement that have positive G-parity. We note that the hermitean 
ordering 



f^af^^ — + (2-15) 



causes the axial-singlet current to vanish, in agreement with its G-parity transformation 
property. Note also that the expression (p.l4|) is of order l/N^j and thus not of leading order 



in the 1/Nc expansion. 

Thus the studies in the simple Skyrme model indicate that an hermitean ordering resolves 
the ambiguities, which occur in the transition from the classical to the quantum level, in 
agreement with the constraints imposed by the G-parity symmetry. 

3. Including Vector Mesons 

It is well known that the inclusion of vector mesons improves many predictions of the 
simple Skyrme model. The too small value of qa is slightly increased but the problem 
remains. Therefore one is naturally curious about possible 0{1) corrections to the leading 
0{Nc) expression for (see table |2]l|). In this section, the G-parity of the currents in 
the model with vector mesons is studied. It is shown that if an arbitrary operator ordering 
were to be allowed then 0{1) corrections to both and the isovector magnetic moment, nv 
would appear to exist. However, such contributions are seen to be G-parity violating and 
hence must be eliminated by an appropriate operator ordering. 

In the presence of vector mesons the situation is more involved not only because the static 
energy functional, Ec\ contains these additional fields but also because field components, 
which vanish classically, get excited by the collective rotation ( p.4| ). At this point we will not 
present the detailed structure of the Lagrangian we are considering. Its form has already been 
proposed a decade ago Also the static soliton configuration has been constructed 



and the fields induced by the collective rotation have been computed |]22[ . These allowed for a 
reasonable description of static baryon properties even in the three flavor model Here we 
rather wish to discuss the G-parity properties of the currents. Besides the pseudovector p^, 
which is defined after eqn ( |2.13| ), we need the isovector- vector i;^ = ^^9^^ — <9^^^^ Obviously 
Vfj, is even under G-parity. The vector meson fields are most conveniently parametrized in 
terms of 

i 

Ri, = ^,, + Pf^- —Vf,. (3.1) 

The coupling constant g is determined from the decay width of the process p — ^ tttt. As the 
isoscalar-vector field u is odd under G-parity no definite G-parity quantum number can be 
attached to i?^. The isovector- vector field possesses the same quantum numbers as v^. 
In terms of these quantities (and the field tensor for the vector meson fields) the symmetry 
currents of the model under consideration are displayed in eqn (2.12) of ref. ||2^. There are 
many terms so we will explicitly display here just two representative ones for the axial vector 
current: 
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Table 3.1: The behavior of the vector meson radial functions under F —F and the G- 
parity quantum number of these fields. 







Hr) ^i{r) ^2{r) 


F ^ -F 
G-parity 


+ 
+ 


+ + 
+ + 



+ pV^" - p'^R^p'') ^] I + • • • ■ (3-2) 
Here my denotes the vector meson mass while the coupling constant 71 can be determined 



from anomalous decays like lu tttttt I^T]. Since is given by the matrix element of 
/ d^r 2Al between proton states with spin projection +| we only need to consider the spatial 
components . In order to discuss the G-parity properties of A'}, we need both the static 
as well as the excited vector meson fields. The static hedgehog ansdtze read 

1 G(r) 

^oi'^) = TT^ir) and p^(r) = e-ikafk (3.3) 

2g gr 

while all other field components vanish. The equations of motion for the radial functions 
Lj{r) and G{r) as obtained from minimizing the classical energy functional are, in general, 
non-linear inhomogeneous differential equations. The source for uj{r) is of the form Vi(r) 
(c/. table pTip while G{r) ~ eika^k tr(rat>j). Hence uj{r) changes its sign under the reversal 
F — >■ —F but G{r) does not. Thus we find that again the transformation F — >■ —F gives the 



proper G-parity quantum number. Upon the collective isorotation ( p. 41) we have not only 



p"(r) A(t)pi{r)A'' (t), but, as already mentioned, additional fields are induced 

c^. = ^e,,k^,h and po = ^^W^ " [^i(^)^ + ^2(0(^ • ^)r])AKt)- (3.4) 

It is important to remark that the ordering between A{t) and f2 (— > J /a^) is completely 
arbitrary because they are considered to be commuting c-number quantities. The radial 
functions ^i(r), ^2(?") and $(r) solve linear inhomogeneous differential equations which serve 
to extremize the moment of inertia, a^. The classical fields act as inhomogeneous parts in 
these equations. Hence one can deduce the behavior of the induced fields under F —F. 
It turns out that ^i(r) and ^2('") are even while $(r) is odd, in agreement with G-parity 
of these fields. These results are summarized in table It is, of course, not surprising 
that the behavior of the vector mesons under F —F coincides with their G-parity as the 
underlying Lagrangian possesses this symmetry and the vector mesons couple according to 
their isospin and spin quantum numbers to the currents displayed in table 

Now we are completely equipped to discuss possible ordering ambiguities for gA and their 
eventual resolution by demanding the proper G-parity behavior. We will see that up to next- 
to leading order in the 1/Nc expansion no contradiction to the G-parity symmetry appears 
as long as we consider all fields to be classical objects. Due to its isoscalar character the uj 
meson does not contribute in the first term in eqn (|3.2|) , which is "non-anomalous" {i.e. no 
e-symbol). It is thus obvious that this term is odd under G-parity just as the isovector-axial 
vector current is supposed to be. For the anomalous term the situation is somewhat more 
involved. One can convince oneself that for this term the proper G-parity is obtained only 
when R is odd, i. e. only the uj meson is allowed. Since we are taking p = "i, one of the 
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three remaining indices has to be associated with the time coordinate. Let us first assume 
that this is attached to the uj field, which is classical. The p's are classical as well because 
only their spatial components contribute. Thus the total contribution from the anomalous 
terms to A1 is classical and no ordering ambiguities appear when the time like coordinate is 
assigned to the uj meson. This, of course, is nothing but the leading order part of A" in the 
1/Nc expansion. It contributes to the radial functions A^{r) and A^{r) in table pTT| and has 
been studied in detail [^, What happens when the uj field is spatial? In that case the 



contribution of the anomalous term to contains two powers of since one of the two p's 
contains a time derivative. Hence this term is of next-to-next-to leading order in the 1/Nc 
expansion and is commonly neglected. 

A next-to leading order term can be obtained only when the isovector part of is 
considered and the angular velocity is treated as an operator in the collective space. Let us 
see how this comes about. A straightforward calculation assuming the apparently natural 
ordering indicated in (|3.4| ) leads to 

= ^ 71 / [rF' (G + ^i) sinF - (6 + 6) sin^F] ^tr ( y^^s^^) ^ (3.5) 

where the matrix element with respect to the nucleon state is indicated. Classically the trace 
in (|3.5|) vanishes identically. However, if one replaced Vl^ by Jzio? this would no longer by 
the case. Instead one would find by using ( p.8| ) and ( |2.11|) 




However, this result cannot be taken seriously since the integrand in eqn (|3.5|) has the incor- 
rect transformation property when going from the nucleon to the anti-nucleon. According 



to table ^]T] the integrand should be odd under F —F but it is even as can be observed 
from table pTT| . Physically ^ implies that the decay amplitude for the /5-decay of the 
anti-neutron is different in magnitude from that of the neutron. This, of course, contradicts 
our present understanding of physics. Hence we conclude that the collective operators have 
to be ordered such that ( |3.5|) vanishes identically. 

Here we have just considered two terms of the complete axial vector current in the vector 



meson model. The above arguments, however, apply as well to all terms in (2.12) of ref.[23 
Furthermore one may conduct the same studies for the isovector-vector current in this model. 
The spatial components of this current give the isovector part of the magnetic moment /ly . As 
expected, the leading order in the 'VjNc expansion is free of ambiguities and agrees with the 
G-parity constraints. At the next-to leading order ambiguities for the collective operators 
occur. Again the improper behavior of the radial part of the matrix element under F —>■ — F 
forces one to arrange these operators so that the 1/Nc correction to fiy vanishes. 

The lesson we learn from studying the symmetry currents in the vector meson model is 
that matrix elements which are different from zero only because of the commutation relations 
( |2.8D but vanish classically are likely to violate the G-parity symmetry. In that event those 
matrix elements should be discarded. 

4- Chiral Quark Solitons 

The results found in the vector meson model make us suspicious about the recently 
discovered l/Nc corrections to qa and fiy since they depend on operator ordering. 

Especially, it has been noticed that both the vector meson addition and the quark addition 
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to the non-linear a model similarly describe short distance corrections to the pion cloud. 



This is true both for the problems of the neutron-proton mass difference and the matrix 



element of the axial singlet current ||2^ . 

Unfortunately the effect of the reversal F —>■ —F is not as clear in quark soliton models 
as it is in purely mesonic models. We therefore divide the present section into two parts. 
In the first one we set up the chiral quark model for both the nucleon and the anti-nucleon 
sectors. We derive a prescription to transform the quark spinors between these two sectors 
(c/. (|4.9|) ). In the second part this is applied to the investigation of the behavior of the 1/Nc 
corrections to qa and /iy under F —F. Technical details of these studies are presented in 
appendix A. 



The chiral quark model |2g, , the simplest chirally symmetric model to contain quark 
solitons, is defined as the sum of the non-linear cx-model Lagrangian with a pion mass term 



Cnia = ytr (d^Ud^U^) + '-^ti (f/ + f/t - 2 
and a valence quark field \l/vai in the background of the chiral field U 



ml 

TV J '> 



(4.1) 



(4.2) 



Here the mass m = g^f'Tr represents a convenient parametrization of the coupling g^^ between 
the valence quark and the chiral field. In this model the eventual effects of sea quarks are 
assumed to be represented by the kinetic term in (4J.). It is also obvious that the model is 
formulated locally. 

Again we employ the hedgehog ansatz ( p.2|) for the pseudoscalar fields. For this static 



(4.3) 



configuration the Dirac equation becomes an eigenvalue problem 

which defines the Dirac Hamiltonian 

h{F) = OL ■ p + f3m cosF(r) + i'^^T ■ r sinF(r) 



(4.4) 



as a function of the chiral angle, F{r). This Dirac Hamiltonian commutes with the grand 
spin 



G = l 



T 

2 



(T 

2" 



(4.5) 



and the parity (H) operators. Thus the eigenstates are classified by G^'"^*- where Hint refers to 
the intrinsic parity, which is defined via H = Hint(— 1)*^. The construction of these eigenstates 
is described in appendix A. 

The eigenvalue evai turns out to be a functional of F, as does the total classical energy 



E,, = + sgn{B)Nce 



val 



(4.6) 



with Eni, = 27r/2 / dr (r^F'^ + 2sm''F + 2ml{l - cosF)) . The sign of the baryon number 
has been included in order to accommodate the hole interpretation of the Dirac theory. In a 
moment we will see that this relative sign is also obtained by requiring that the nucleon and 
anti-nucleon possess equal masses. 
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The soliton is computed by solving the Euler-Lagrange equation for the chiral angle 

2 sin2^ N tji f dVl 

F" = —F' + —— + m^sinF - sgn(5)^^ / —^l,/3 (sinF - i^,r ■ rcosF) ^.^i (4.7) 

self-consistently. For the boundary condition -F(O) = — vr the bound valence quark state of 
the self-consistent solution is found in the C^'"' = 0''" channel. We denote the corresponding 
eigenvalue and eigenfunction by e^j and ^I/^i, respectively. On the other hand, when assuming 
the boundary condition F{0) = vr, a strongly bound quark is only obtained in the 0~ channel 
with eigenvalue e~^i and eigenfunction This, of course, reflects nothing but the fact 

that the parities of the nucleon and the anti-nucleon are opposite. In fact the self-consistent 
solutions for the nucleon {F{0) = — tt) and anti-nucleon {F{0) = vr) are distinguished by 
different overall signs of F(r) and the eigenvalues of the bound quarks, i.e. e~^i = — e^j. 
Thus, due to the inclusion of the factor sgn(i?) in ( ^^ , the classical masses of the baryon 
and anti-baryon are identical. 

Actually, the self-consistent solution for the anti-baryon can easily be obtained from the 
one for the baryon by noting that 

h{-F) = -J^h{F)J with J = J^ = if3j5- (4.8) 

The transformation J' commutes with G but has negative parity. This implies for the 
eigenvalues and eigenstates of the Dirac Hamiltonian 

ef^^^-ef and \fj.,G^) = J\fx,G^), (4.9) 

where fi labels the particular eigenstate. Details of the transformation ( [4.9| ) in terms of the 
radial parts of the quark wave-functions are given in appendix A. We will make extensive 
use of this transformation when comparing the currents for the nucleon and anti-nucleon. 
However, we first have to perform the projection of the chiral quark soliton onto states with 
good spin and isospin. For the non-linear a model this is straightforwardly achieved by 
substituting the rotating hedgehog ( P^ ) into the defining equation (^.1|). This yields the 
mesonic part of the moment of inertia 



vr 



=- 



J drrhm^F (4.10) 



which obviously is identical for the baryon and the anti-baryon. In order to compute the 
quarks' contribution to we employ the cranking method [^]. Then the quark spinors are 
rotating in isospace, which adds the Coriolis term (1/2)t ■ O as a perturbation to the Dirac 
Hamiltonian ([4.4]). The energy eigenvalue acquires a change in second order perturbation. 



which (due to isospin invariance) may be written as evai Cvai + (l/2)a;qO . Here 



q 

= (4.11) 



denotes the quarks' contribution to the moment of inertia of the nucleon. As [J , Tj] = one 
finds that a^F] changes sign under the transformation ( |4.9| ), aq[— F] = — aq[F]. Taking into 
account that ( [4.11[ ) represents a perturbation to valence quark energy, which enters the total 



energy with the overall factor sgn(i?), the total moment of inertia 

= al + sgn{B)al (4.12) 
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is seen to be invariant under F — >■ —F. This guarantees equal masses for the baryons and 
anti-baryons in the presence of rotational corrections. 

In the context of the cranking method the induced components of the valence quarks are 
computed in a perturbation expansion in powers of O. This will serve to make the discussion 
of the particle conjugation properties of their contribution to the currents more transparent 
than in a variational approach [^. In order to be consistent with the expansion for the 
energy we require the first order expression for the cranked wave-function 

*e„.. = Ait)U., + W W^l^ I . (4. 13) 

The isospin matrices r carry unit grand spin and are positive under parity. Thus G 1^ 
for F(0) = ±7r. 

We now turn to the discussion of the currents in the chiral quark model. These are again 
sums of mesonic and quark parts. As they are the Noether currents of the underlying theory 
the overall factor sgn(i?) is carried along for the quark part 

y^^K) = \f> {I (iV,i^ T i^V.i) ] + Sgn(5)^„,,k7;.(75) Y^crank. (4.14) 

The interesting part is, of course, the one which is due to the quarks, especially the spatial 
components (Q^ = ai{ai) Ta/2) 

A'c4'L„kS!'*=r..k = NcD^¥„,Q\^„, (4.15) 

^ /x^val ^val Z Cval 

In this expression we have temporarily adopted the ordering between Dab and O suggested 
by the form of the perturbed wave-function ( 4.13 ). In the framework of the chiral quark 



model the 1/Nc corrections to qa and /xy, which were recently discussed in various similar 
models H, H, 0, |[, correspond to the replacement f2 — ^ J/a^ in the ordering of ( [4.15| ). 

To obtain qa we calculate the integral / d^rA^. The contribution from the next to leading 
order involves the sum over the grand spin projection (in which the energy eigenvalues are 
degenerate) 

^(val,±|a3r,|/i,lM,T)(/x,lM,T|r,|val,±) = '-S^L^.e^by (4.16) 

The signs appearing in the state vectors label the intrinsic parity. The upper (lower) sign 
refers to the (anti) nucleon sector. We refer to appendix A for the actual computation of 
these matrix elements and the definition of the integrals on the RHS. Since these matrix 
elements are proportional to the anti-symmetric e-tensor the appearance of the commutator 



[Dabi Jj] is obvious (see (|2.8|) ) when replacing the angular velocity by the spin operator. 
Taking account of this prescription one finds for the axial charge of either the nucleon or the 
anti-nucleon 
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where the 1/Nc correction, g^^ , is represented by the term proportional to We should 

remark that ( [4.17|) corresponds to a spin up proton or a spin up anti-neutron which each give 
the same matrix element of D33. Note that the overall sign differs in the two cases. Similar 
calculations give for the isovector part of the magnetic moment 

^^ = ^f+^i^) = ^-IflMj, j drrhin'F (4.18) 
-sgn(5)— fM^ (val| r-a,a3 r3|val) + — ^ ^ ^'^/^^l ^ ' ' ^ 

Here the e-tensor, which again causes the appearance of the commutator ( p.8| ) in the next-to 
leading term, originates from the definition of the isovector part of the magnetic moment 
/iy = l/2((r X V^). 

We have computed these matrix elements in the chiral quark model numerically. It turns 
out the next-to leading order is only about 1/10 of the leading order. Thus these terms are 
more strongly suppressed than in the NJL model. This is apparently due to the fact that the 
chiral angle in the chiral quark model has quite a large extension yielding a larger moment 
of inertia. 

Now we can investigate these expressions with respect to their behavior under particle 
conjugation. This is made simple by noting that we just have to apply the transformation 
( ^^ . Since Jai = —aiJ and Jai = OiJ it is obvious (noting the factor sgn(i?)) that the 
leading order terms transform properly in agreement with table ^]l|, i. e 



(val, +10-3X3 1 val, +) — > (val, -|o-3r3|val, -), (4.19) 
(val, +1 [r • Q;,o-3]r3|val, +) — > -(val, -| [r ■ a, 0-3] rgjval, -), (4.20) 

Since the transformation J does not affect the isospin structure of the matrix elements it 
is obvious that the matrix elements involved in the terms describing the 1/A^c corrections 
transform in the same way 

(/i, -|ri|val, +)(val, +|(T3r2|/i, -) — > (/i, +|ri|val, -)(val, -|(J3r2|yU, +) (4.21) 

p ^ p 

(/X, +|ri|val, -)(val, +| \r ■ a, 0-3] r2|/i, -) — > -(/i, +|ri|val, -)(val, -| [r • a, 0-3] r2|yU, +). 

In appendix B we present an additional argument, which is based on the completeness of 
the eigenstates of the Dirac Hamiltonian, for the transformation relation ( |4.21| ). Now the 
G-parity violation of the ^jNc corrections becomes clear: The energy denominator in these 
terms changes its sign when going from the nucleon to the anti-nucleon. Thus the 'VjNc 
corrections transform oppositely to the leading order contribution, i.e. incorrectly. We 
therefore conclude that the ordering ambiguities contained in the formulation of the quark 
currents ( [4.15|) have to be resolved in such a way that the 1/A^c corrections to qa and [ly 



vanish: 

c/S^ = and [i^^ = 0. (4.22) 
This can be achieved by demanding an hermitean ordering prescription 

Dab^j — . ^ {DabJj + JjDab) • (4.23) 
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For this ordering each of the two subleading terms in ( [4.14| ) vanishes in the nucleon subspace 
after performing the sum over the grand spin projection ( |4.16| ). Since 



4 

ekbj {DabJj + JjDab) = -^hekbj Wj + JjJb) = (4.24) 



the required result (|4.22| ) is obtained. 



It is seen that the pattern obtained in the vector meson model is repeated in the chiral 
quark model. As already noted we are considering a local formulation of the chiral quark 
model. Non-local effects, as e.g. the blocking of transitions to negative energy states in 
( [4.13 ), associated with the inclusion of levels other than |val) in the definition (O) may 



cause different expressions for g^^^ and fiy\ We will comment on related studies 0, below. 

To avoid confusion, we should remark that there is no practical problem in the chiral quark 
model with too small gA- In fact, gA comes out slightly too large in this model, although the 
problem can be resolved [^by complicating the model a bit. However, the structure of the 
chiral quark model is, for present purposes, the same as other quark soliton models. 

Let us conclude this section by commenting on the G-parity behavior of the axial singlet 
current in the chiral quark model. This has been noticed some time ago to account for the 
proton spin puzzle fairly well When the isospin part of the generator reduces to a 
unit matrix, the first term on the RHS of ( |4.15| ) is easily observed to vanish. Furthermore no 



ordering ambiguities occur in the second and third terms because the D-matrix is replaced 
by the unit matrix. The radial part of the matrix element between nucleon states of the axial 
singlet current 

/ d'rAl = sgn(5)i + 

— ' sgn(5)— Y Z ('"l^slval) + h. c. , (4.25) 

is invariant under F —F. This result indicates that the induced part of the quark current 
behaves properly under G-parity if there happens to be no ordering ambiguity and the leading 
term in the 1/Nc expansion vanishes. 

5. Summary and Conclusions 

We have formulated the particle conjugation operation, and especially its convenient 
realization as G-parity, in the framework of chiral soliton models. We specifically treated 
the usual Skyrme model, the Skyrme model with vector mesons and, as the simplest example 
containing quark fields, the chiral quark model. The prescription for G conjugation in the 
solitonic sector is simply to reverse the sign of the profile function F{r) —F{r) while 
leaving the collective coordinates A{t) unchanged. 

It was found that requiring the correct particle conjugation properties for currents could 
resolve operator ordering ambiguities which greatly affect the results of a number of interest- 
ing calculations. These ambiguities occur in the transition from the classical to the quantized 
formulation in the space of the collective coordinates. 

In the basic Skyrme model we noted that an hermitean ordering of operators in the 
collective space was required to guarantee the vanishing of a G-parity violating axial current 
form factor. The same ordering also guaranteed the vanishing of the G-parity violating 
matrix element of an axial singlet current of interest in connection with the "proton spin 
puzzle" . 



15 



In the Skyrme model with vector mesons we investigated possible corrections down by 
1/Nc to qa and iiy These turned out to vanish to insure G-parity conservation. Generally 
speaking, our studies indicate that objects which vanish classically due to G-parity should 
also vanish in the quantized formulation. 

In the chiral quark model the situation turned out to be somewhat more involved because 
the quark states are not eigenstates of the particle conjugation. Having set up the formulation 
for the anti-nucleon we were able to show that the next-to leading order 1/Nc corrections to 
qa and /xy, which are non-zero for the ordering suggested by the ansatz for the cranked quark 
fields (|4.13|) , behave oppositely to their leading order counterparts under F —>■ —F. From this 



we concluded that the operator ordering in the space of the collective coordinates should be 
arranged so that these corrections vanish. It should be stressed that this conclusion follows 
immediately after ( [4.15|) without reference to the quark wavefunction phase convention. This 



is because the identities ai{ai)J' = ^Jai{ai) and [J,t] = immediately show that the 
leading and next-to leading order piece of the quark currents in the 1/Nc expansion transform 
oppositely. Since the leading order has the correct G-parity behavior this type of next-to 
leading order corrections should actually be absent. This can be accomplished by adopting 
the hermitean ordering prescription (|4.24 ) when going from the classical to the quantum 



level. 

The results we have obtained in the chiral quark model are based on an ordinary pertur- 
bation expansion which is essentially equivalent to the cranking procedure of ref. [^. Using 
a somewhat different expansion scheme for the Dirac operator including the Coriolis term 
the authors of refs. |^ obtained a somewhat different result^ for the 1/Nc correction to 
the axial charge 

V- ^ . (val|ri|/i)(/i|cr3r2|val) 
-sgn S — ^ 2^ sgn . 5.1 

Although this result deviates from the correction in ( ^.17] ) only slightly numerically it is 
obvious from our discussion that ( |5.1| ) transforms properly under the particle conjugation 
due to the additional factor sgn(e^)[^. It should be remarked that the treatment of refs. 

is not itself free of ordering ambiguities, which for example manifest themselves in 
the violation of the PC AC relation. Furthermore the expression ( |5.1|) cannot be obtained 
from the ordinary perturbation expansion or the cranking procedure of ref. It seems 

to us that a number of open questions remain as to the validity of ( p.l|) . The discussion 
of these goes beyond the scope of this paper nevertheless our considerations on the particle 
conjugation symmetry favor the results of refs. 0, |^ over those of refs. |^. The latter are 
essentially identical to ( [4.17| ). Stated otherwise, the considerations concerning the particle 



conjugation symmetry rule out the results of refs. 0, § while no definite statement can be 
made on those of refs. |^ ^. As the models [^, |^ contain non-local effects they are somewhat 
different from the chiral quark model considered here. 

The considerations presented in this paper seem to exclude the emergence of 1/Nc cor- 
rections to Qa and fiy in the context of the collective quantization of the static soliton 
configuration when performing a "classical" expansion in the angular velocity O in local 
models. This suggests that the incorporation of pion fluctuations off the soliton is important 
in order to obtain non-vanishing corrections. Such a treatment seems quite appealing since 
it has been shown to solve problems PT| related to the fact that the current algebra relations 



cannot be correctly described within the collective quantization of the static soliton p2|. 



"We are grateful to C. Christov for pointing out this difference to us. 
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Appendix A: Explicit Solution to the Chiral Quark Model 

In this appendix we present the explicit form of the eigenfunctions of the Hamiltonians 
h{F) and h{—F) for the nucleon and anti-nucleon, respectively. The Hamiltonian h{F) is 
given in (|4.4| ). Also the expressions for the axial charge, qa, the isovector part of the magnetic 
moment, nv, as well the corresponding "would-be" 1/Nc corrections g^^^ and /iy'* are made 
explicit. 

Although only states with grand spin G = and G = 1 are relevant in the framework of 
the chiral quark model, it is helpful to define general eigenstates of the grand spin operator 

\ljGM) = E^Si3E^£|JMI^^3)4^3)z. (A.l) 

j3S3 mis 



Here Cjj^fj'j'^ denote the SU(2) Clebsch-Gordan coefficients. The subscripts J and / indicate 



-iJM 
' 333,3' 

spinors in coordinate- and iso-space, respectively. 

1. Treatment of the Nucleon 

The quark state which is bound in the chiral background field F{r) resides in the (7 = 
channel and possesses positive parity when the boundary condition F{r = 0) = — vr is chosen. 
In terms of the grand spin states ( [A.lD this valence quark state is parametrized as 



ltgt{r)\Q\m) 

|val,+)= |. (A.2) 

\U{r)\i\m) 

The radial functions gQ{r) and /o('") are obtained by diagonalizing the Dirac Hamiltonian 
( [4.4|) and identifying the state with the energy eigenvalue in the interval —m < e^,! < m as 
the valence quark. Technically the diagonalization is achieved by constraining the system to 
a spherical box of finite radius D and imposing certain boundary conditions at r = D. Here 
we demand that the upper components of the Dirac spinor ( |A.2| ) vanish at r = D. For more 



details on the numerical treatment we refer to ref . ||34|| . 

The matrix elements relevant for the classical parts of the axial charge and the iso-vector 
part of the magnetic moments are easily obtained to be 



[r? 



(val, +|a3r3|val, +) = ^ J drr^ (^g+{rf - 

(val, +1 [r ■ a, ^3] r^jval, +) = jS^b J drr^g+{r)f+{r) (A.3) 

In order to evaluate the perturbated wave-function due to cranking we need to compute 
the overlap matrix element (/i|rj|val, +). The isospin matrix tj carries unit grand spin and 
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positive parity. Thus non-vanishing matrix elements exist only when belongs to the 1" 
channel. These states are parametrized as 



l/i, IM, -) 



/z^7i-(r)|2|lM) 
v/r^(r)|lflM), 



ig^^{r)\Q\lM) 
-/2-(r)|lilM) 



(A.4) 



The corresponding eigenvalues of h{F) are denoted by e~. Furthermore M = —1,0, 1 refers 
to the projection of the grand spin G = 1. Again we refer the interested reader to ref. [p4| 
on details on the numerical evaluation of the radial functions gt^-, f2^- The boundary 
conditions on these radial functions are such that the upper component of the spinor 



vanishes at r = D in order to avoid finite size isospin violations ||35|| . 

It is then straightforward to compute the matrix elements relevant for the axial charge 
and iso-vector part of the magnetic moments, which are discussed in section 4: 



(/i, IM, -|rj|val|,+) 
(val|, +|a3rb|/i, IM, -) 

(val, +1 [r ■ a, as] Tfel/x, IM, -) 
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IM 



^12 r I 21 r 



^12r I ^21r 

Th 0-iM + Tf^ OiM 



(A.5) 



The coefficients 5*^, and are given as integrals over the radial functions defined in 
eqns (|A.3| ) and (|A.5|) : 



Si. 



R1. 



drr'^ 
drr^ 



(A.6) 



The summation over the grand spin projection yields 



^(val, +|(T3rf,|/i, IM, IM, -|r,-|val, +) = -S'+L+esbj 



(A.7) 



M 



M 



^(val, +1 [r ■ a, 0-3] n\fi, IM, IM, -|r,-|val, +) = ^S^R^e-^hj 



2. Treatment of the Anti-Nucleon 

For the boundary condition F{r = 0) = vr the parity of the bound valence quark gets 
reversed. Thus this state is member of the 0~ channel: 



|val,-) 



■ig^{r)\l\m) 
/o'WlOiOO) 



(A. 



The corresponding eigenvalue of the Dirac Hamiltonian ([4.4|) is labeled e^ai- As the trans- 
formation JT" in eqn ( [4.8|) exchanges upper and lower components of a Dirac spinor, we now 
demand the lower components of ( |A.8| ) to vanish a,t r = D. This makes possible the compar- 
ison with ( [A.2| ) even for finite D. 
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The classical parts of gA and jjy are given by 

(val, -[asrslval, -) = J drr^ (^f^ {rf - ^g^ {rf 



(val, -| [r ■ a, 0-3] rb|val, -) 



^36 / drr^fo {r)gQ (r) 



(A.9) 



Under the cranking perturbation the valence quark state ( |A.8| ) couples to the states in the 
l"*" channel 



IM, +) 



/ zgUrMlM) 



'igUr)miM) 



\-ft,ir)\2llM)J V/2+(r)|0ilM) 



(A.IO) 



Their eigenvalues are denoted by e+. Again the lower components are constrained to be zero 
at r = D. The isospin structure of the matrix elements ( [A. 51 ) remains unaltered when we 
transfer them to the anti-nucleon. However, the coefficients S^, and i?^ are replaced by 



s: 



R~ 



— j drr"^ 
drr^ 



9o {r)g^^{r) + /g {r)f+^{r) 
3/o-(r)/2+ (r) - goir) {gt.ir) - V2gt^{r 



(A.ll) 



V2g 



When the summation over the grand spin projection is performed these integrals have to be 
substituted for their analogous expressions in eqn (|A.7| ). 

Using the explicit representations (|A.2| , A. 4 , A.8| ) and (|A.10|) the transformation (|4.9| ) 
corresponds to 



9ti.r) -foir), foir) go{r), 

91 (r) ft ir) , /f (r ) ^ -c/+ (r ) , 

92 (r) - ft (r) , ft ij-)^ gt (r) . 



(A.12) 



Numerically we have verified that the self-consistent solution reverses its sign {F —F) 
when we alter the boundary condition from -F(O) = — vr to -F(O) = tt. Furthermore the 
computations confirm the transformations e^j — Cvai ^'^d ^/I ~^ ^'^^ eigenvalues of 
h{F). The transformation properties for the wave-functions ( |A.12| ) are regained up to an 
overall, irrelevant phase. This phase may, of course, be different in the grand spin zero and 
one channels. 

Appendix B: Completeness Argument 



Using the results from appendix A the transformation relations ( [4.19| ) and ( |4.21| ) are 
easily verified. In terms of the integrals over the wave-functions these transformations are 
expressed as 



C+T + 



and StRt 



(no sum over /i) 



(B.l) 



The validity of the relations ( [4. 191 ) and ( [4. 21] ) is also nicely confirmed by making use of the 
completeness property of the eigenstates of the Dirac Hamiltonian. From that one can deduce 

1 i 
(val, ijasTslval, ±) = - ^t^t ±1 ['^ ' ^3] Valval, ±) = - S^^t- (^-2) 
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This indicates once again that the matrix elements involved in the subleading order of the 
1/Nc expansion transform identically to those of the leading order when the sign of the chiral 
angle is reversed. We should mention that we have confirmed the identities ( [B.2| ) numerically 
for both boundary conditions -F(O) = ±7r. 
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